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ABSTRACT 

How do peculiar velocities affect observed voids? To answer this question we use 
the VIDE toolkit to identify voids in mock galaxy populations embedded within an 
N-body simulation both with and without peculiar velocities included. We compare 
the resulting void populations to assess the impact on void properties. We find that 
void abundances and spherically-averaged radial density profiles are mildly affected by 
peculiar velocities. However, peculiar velocities can distort by up to 10% the shapes for 
a particular subset of voids depending on the void size and density contrast, which can 
lead to increased variance in Alcock-Paczyhski test. We offer guidelines for performing 
optimal cuts on the void catalogue to reduce this variance by removing the most 
severely affected voids while preserving the unaffected ones. In addition, since this 
shape distortion is largely limited to the line of sight, we show that the void radii are 
only affected at the ^ 10% level and the macrocenter positions at the ^ 20% (even 
before performing cuts), meaning that cosmological probes based on the Integrated 
Sachs-Wolfe and gravitational leasing are not severely impacted by peculiar velocities. 

Key words: stacked cosmic voids—peculiar velocity—simulations—large-scale struc¬ 
ture of universe. 


1 INTRODUCTION 

In the era of the large scale surveys, extracting cosmological 
information from the large scale structure of the Universe is 
a major challenge. The distribution of galaxies in the Uni¬ 
verse follows a cluster-sheet-void-filament pattern called the 
cosmic web (Bond et al. 1996). Recently, the underdense re¬ 
gions of the cosmic web (discovered in 1978, see Gregory Sz 
Thompson (1978)) captured the interest of cosmologists, as 
they open up new avenues of constraining models of the evo¬ 
lution of the Universe. Cosmic voids promise to be dynam¬ 
ically simpler than overdense structures (van de Weygaert 
& van Kampen 1993; Hamaus et al. 2014) and fill most of 
the Universe (Fairall 1998; Hoyle & Vogeley 2002; Bos et al. 
2012): they are potentially powerful tools to constrain cos¬ 
mological models. 

Voids are being used for a wide range of applications, 
spanning from their use as standard spheres through the 
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Alcock-Paczyhski test (Sutter et al. 2012a; Sutter et al. 
2014c), to the study of the late-time Integrated Sachs-Wolfe 
effect in f{R) models (Cai et al. 2014a), or to constrain 
models of modified gravity (Li 2011; Spolyar et al. 2013; 
Clampitt et al. 2013; Zivick et al. 2014; Cai et al. 2014b), 
models of coupled dark energy (Sutter et al. 2015) and cos¬ 
mological models of massive or massless neutrinos (Massara 
et al. 2015). 

In this context, it is crucial to properly model the im¬ 
pact of systematic effects, which might undermine the qual¬ 
ity of cosmological constraints from cosmic voids. Voids 
probe a different dynamical regime of structure formation 
than high density regions, they show a quasi-linear be¬ 
haviour (van de Weygaert V van Kampen 1993; Hamaus 
et al. 2014) and are surrounded by mildly non-linear walls 
(Paz et al. 2013). For this reason, they present different sys- 
tematics than cosmological tools from high-density regions. 

The peculiar velocities of galaxies are the main sys¬ 
tematic when using voids for cosmology (as discussed by 
e. g. Ryden (1995) and Sutter et al. (2014c)). The position 


© 2015 RAS 



2 A. Pisani et al. 


of each galaxy is measured in redshift space, thus the pe¬ 
culiar velocity of the galaxy changes the galaxy position by 
adding a component to the distortion due to the expansion 
of the Universe — hence affecting our measures of the cos¬ 
mic web. It is therefore important to assess the impact of 
such velocities on voids when aiming to a cleaner and opti¬ 
mal extraction of the cosmological signal. Recently Lavaux 
& Wandelt (2012) studied the statistical impact of velocities 
on a dark matter simulation, and found a uniform bias in 
the Alcock-Paczyhski measurement when peculiar velocities 
were included. This was confirmed by Sutter et al. (2014c) 
for mock void populations, but those works did not examine 
in detail the source of that bias. 

Other works focused on avoiding such bias, for example 
Pisani et al. (2014) presented a method to extract informa¬ 
tion from voids in a model-independent way through the 
reconstruction of voids density profiles in real space, that 
promises to improve the modelling of systematics for future 
surveys. Also, Hamaus et al. (2014) introduced void auto¬ 
correlations as a promising statistic with reduced systemat¬ 
ics impact. Nevertheless in most cosmological uses of voids, 
we lack of a detailed analysis of the velocity effects. While a 
full treatment of such effects on voids could lead to a reduc¬ 
tion of systematics, its implementation is non-trivial (Paz 
et al. 2013). 

Thus, in this work, we focus on the effect of velocities on 
void statistics and properties using mock galaxy catalogues, 
to make contact with voids measured with current surveys 
and used to constrain cosmology (as in Sutter et al. (2014c)). 
The analysis of such effects is instrumental to a greater un¬ 
derstanding of void systematics and avenues for correcting 
for their impacts. Rather than only examining a single void 
in isolation we study an ensemble of voids evolved in a full 
simulation, analysing their properties both on a one-to-one 
basis and on average. 

Works based on cosmological analysis of voids from 
real data (e. g. ISW (Cai et al. 2014a), leasing (Melchior 
et al. 2014), Alcock-Paczyhski test (Sutter et al. 2014c), 
abundances (Pisani et al. 2015), which could be affected 
by peculiar velocities, can use our results as a measure of 
the dynamical effects affecting the observed voids, and con¬ 
sider the guidelines emerging from this work to exclude the 
most affected voids from the analysis. Using such guidelines, 
the signal-to-noise ratio for cosmological measurements such 
as the Alcock-Paczyhski test (Alcock & Paczynski 1979) 
promises to be enhanced. 

The paper is organised as following: in Section 2 we 
present the void finder and simulation employed for this 
work, as well as the matching algorithm and methodology 
used to compare the two void catalogues (with and without 
velocities). In Section 3 we present the results: we discuss 
which voids survive after the application of peculiar veloci¬ 
ties and analyse their features. In Section 4 we focus on one- 
to-one void comparisons to measure the impact of peculiar 
velocities on void properties, such as relative ellipticities, 
radii, and macrocenter positions. We then analyse the effect 
on void stacks, widely used for the Alcock-Paczyhski test 
and for applications such as the integrated Sachs-Wolfe ef¬ 
fect (ISW) or leasing; and on abundances. Finally, in Section 
5 we conclude with a discussion on how to reduce the impact 
of systematics related to peculiar velocities, we comment the 
results and their applicability for current and future surveys. 


2 ANALYSIS METHOD 
2.1 Void finder 

The void finder used for this work is VIDE, an improved ver¬ 
sion of ZOBOV (Neyrinck 2008). The VIDE toolkit (Sutter et al. 
2014a) finds cells in a distribution of tracers by means of the 
Voronoi tessellation and defines basins joining cells through 
the watershed transform. Following ZOBOV’s methodology, 
VIDE does not merge zones if the minimum density along 
the ridge between zones is higher than 0.2 times the mean 
density of the simulation. As illustrated by Neyrinck (2008), 
this criterion prevents voids from growing into haloes. VIDE 
provides the full void hierarchy, but for this work we will 
consider only parent voids with central densities less than 
0 . 2 p (see Sutter et al. (2014a) for details). Other choices may 
also be considered (see, for example, Nadathur & Hotchkiss 
(2014)); according to the discussion in Lavaux & Wandelt 
( 2012 ) and Sutter et al. (2012a), we focus on this choice to 
measure the impact of velocities in the framework of recent 
applications of the Alcock-Paczyhski test. 

With the convention used by VIDE, the total volume V 
of the void is the volume of all the Voronoi cells composing 
the void. The effective radius Reff of the void can thus be 
defined through the volume as the radius of a sphere with 
volume V: 

/ 3 \ 1/3 

«•« = (■) 

We also define the void macrocenter: it is the volume- 
weighted center of the void Voronoi cells. We remind that 
VIDE naturally excludes voids with effective radius below the 
mean particle separation. This allows to set a limit for voids 
of too small size, affected by shot noise. Finally, we recall 
that VIDE computes the shape of a void by taking the Np 
void member particles and constructing the inertia tensor: 

Np 

M,, = Y^ixl + zf) ( 2 ) 

i = l 

Np 

^xy — ^ ^ Xil/i 

i=l 

where Xi, yi, and Zi are the coordinates of the particle 
i relative to the void macrocenter. Defining Ji and J 3 as 
the smallest and largest eigenvalues of the inertia tensor, we 
obtain the definition of the ellipticity of the void: 



The macrocenter, the ellipticity and the effective radii 
of voids allow to analyse the properties of voids. 

2.2 Simulation and HOD details 

The simulation we use in this work is a 1 h“^Gpc box size 
dark matter N-body simulation, for which accuracy and er¬ 
ror behaviour have been improved using the 2HOT code 
(Warren 2013) for cosmological volumes. It contains 1024^ 
particles and has a particle resolution of 7.36 x 
We used 2LPTIC (Crocce et al. 2006) and CLASS (Bias et al. 
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2011) to generate initial conditions. The 2H0T code opera¬ 
tion scales as NlogN in the number of particles. More details 
on the simulation can be found in Sutter et al. (2014b). 

To obtain two mock catalogues, we apply the Rocks tar 
halo finder and use it as an input for an Halo Occupation 
Distribution model (Tinker et al. 2006; Zheng et al. 2007). 
The model assigns to each dark matter halo of mass M a 
central galaxy and satellite galaxies, the mean number of 
central galaxies and satellites is described by: 


{iVcen(M)) = 


1 _l_ erf r ~ logMiiin 


\ CTlogM 

(iVsat(M)) = (iVeen(M))( ^^,-^° 


( 4 ) 


where we have aiogM, Mmin, Mo, Mi and a as free parame¬ 
ters which are set to match the properties of a given galaxy 
population (see Sutter et al. (2014b)). 

We produce two galaxy catalogues, to mimic a high- 
and a low-resolution galaxy sample (with a tracer density 
of 2 X 10 ^ per cubic h ^Mpc for the HighRes sample and 
3 X 10“^ per cubic h“^Mpc for the LowRes sample). We 
assign the halo velocity for each galaxy, and additionally 
give a random peculiar velocity to each satellite galaxy by 
drawing from a Maxwellian distribution with mean equal to 
the dark matter halo velocity dispersion. We then use these 
galaxy samples to study the effects of peculiar velocities on 
the finding of voids in a realistic situation by comparing 
mocks with and without the presence of velocities in the 
high and low density case. 


2.3 The matching algorithm 

To analyse the effects of peculiar velocities on the void de¬ 
tection, we use the matching algorithm provided with the 
VIDE toolkit, able to compare two voids catalogues in the 
most efficient way for our purposes. 

A simple but crucial point for our analysis is the choice 
of the catalogue with peculiar velocities as a base catalogue 
for the comparison. This choice is instrumental, since it con¬ 
siders the correct perspective for the void finding: when we 
find voids in real surveys we observe them with peculiar ve¬ 
locities. Thus any study of peculiar velocities has to use the 
information we have as a starting point, so that the results 
of the analysis can be applied to a real galaxy survey, where 
we only have voids with peculiar velocities. 

For each void in the peculiar velocity catalogue, the 
matching algorithm selects all possible match candidates 
with centres lying in the Voronoi volume of the void. 
The matching method uses the unique particle IDs to find 
matches; it takes as the best match the candidate void with 
the highest number of shared particles. 

We use the matching methodology to perform two re¬ 
lated analyses: first we apply the matching algorithm to 
check which voids found in the peculiar velocity mock cor¬ 
respond to voids in the mock without velocities (see Section 
3). The unmatched voids are highly affected by peculiar ve¬ 
locities, at the point that they do not clearly correspond to 


any real-space void counterpart^and may be created by the 
effects of peculiar velocities. 

With the aim of characterizing the unmatched voids, 
we analyse their properties (such as ellipticity, density con¬ 
trast or radius). Finally we disregard voids that are in the 
catalogue without peculiar velocities but not found in the 
catalogue with peculiar velocities, because in any way we 
will never be able to detect them — the fact that we do not 
find them in the peculiar velocity mock means that velocities 
have erased them. Then we follow this analysis with a sec¬ 
ond study based on the relative properties of the matched 
voids (see Section 4), to measure the impact of velocities 
on a one-to-one basis. To assess the impact of velocities we 
thus consider the unmatched voids, as well as the properties 
of the matched voids, that is voids we measure in redshift 
space corresponding to voids without peculiar velocities. 

The considered approaches are targeted to provide 
guidelines on the effect of peculiar velocities to all applica¬ 
tions using voids from real data. As discussed above, for such 
applications, to avoid Poisson noise effects, voids with radii 
below the mean particle separation are already excluded, 
namely lower than 8 h“^Mpc for the HighRes sample and 
15 h“^Mpc for the LowRes density sample. 


2.4 Example comparison 

To obtain a high quality one-to-one void comparison we use 
the matching routine of VIDE described above. Nevertheless, 
it is very beneficial to obtain a good visualisation of voids, 
to check the behaviour of the matching algorithm and to 
understand the behaviour of voids when peculiar velocities 
are added. 

Since peculiar velocities will impact the detailed shape 
of voids, we represent them visually using the Voronoi cells 
constituting the void. Each cell is represented as a sphere, 
its area is related to the volume of the cell, and the center 
of the sphere is the particle position. The algorithm able to 
represent voids following this idea is part of the public void 
finder VIDE (Sutter et al. 2014a)^. 

Although the Voronoi cells are not spherical, this consti¬ 
tutes an approximation that allows to observe the shape of 
voids in a particularly effective way. Figure 1 shows the abil¬ 
ity of the matching algorithm: it gives a visual impression 
of the quality of the matching on individual voids, allowing 
to see the effect of peculiar velocities on a one-to-one basis. 
The representation shows the amount of shape variation of 
voids, thus serving as a guide for the analysis of the peculiar 
velocities effect: we observe the slice and the matching voids 
from the two different catalogues. 

To illustrate the effect of velocities along the line of 
sight, we show in Figure 2 examples for the effect of velocities 
for a few voids. In this case, the line of sight is towards the 
left, i. e. we represent the voids on the y — z plane. In Figures 
1 and 2 we qualitatively observe minor perturbations on the 
detailed void perimeter, for these particular example voids 


^ With the void definition we are using. Of course, with a different 
definition of voids, such as the one proposed by Lavaux & Wandelt 
(2010), that is based on a Lagrangian orbit reconstruction, the 
situation might change. 

^ At http://bitbucket.org/cosmicvoids/vide_pubIic. 
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^ 54 fe * Mpc (HOD HighRes) ^ 54 hr^ Mpc (HOD HighRes) ^ 54 fe ' Mpc (HOD HighRes) 



-150-100 -50 0 50 100 150 -150-100 -50 0 50 100 150 -150-100 -50 0 50 100 150 


Figure 1. Visual impression of the matching procedure for void catalogues. Left frame shows a slice of the density field without peculiar 
velocities. Central frame shows the same slice but with peculiar velocities and the void found in the slice (red). We represent the Voronoi 
cells constituting the void as small spheres, each with area related to the volume of the cell and center defined by the macrocenter of the 
cell. A visual comparison of the density field between the two panels illustrates the effect of peculiar velocities: as expected, structures 
are slightly enhanced. Right panel shows the void found in the density field without peculiar velocities (blue). The representation of voids 
is on the x — y plane, units are h“^Mpc, and the slice is along the line of sight. Top row represents a void from the HighRes sample, 
bottom row a void from the LowRes sample. 


-)• 63 V^Mpc (HOD HighRes) 



-100 -50 0 50 100 


-> 78 h-^ Mpc (HOD HighRes) -> 87 Mpc (HOD LowRes) 



-100 -50 0 50 100 -150-100 -50 0 50 100 150 


Figure 2. Slice of the density field in the y — z plane, the line of sight is thus towards the left, qualitatively showing the effect of peculiar 
velocities on voids; the slice is with peculiar velocities. According to the color code in Figure 1, we show the void found in the density 
field with peculiar velocities in red, and the void found in the density field without peculiar velocities in blue. Left and central panel 
represent voids from the HighRes sample, right represents a void from the LowRes sample. We see that for these particular voids peculiar 
velocities only affect the detailed outer structure of the void, with minor effect on global properties such as radius and position. 
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the effect of peculiar velocities on global properties such as 
radius and position is negligible. 


3 RESULTS: VOID SURVIVAL RATE 

As discussed in the previous section, we used the catalogue 
with peculiar velocities as a base catalogue for the compari¬ 
son. In order to extract cosmological information, we would 
like to know which voids are the most affected by peculiar 
velocities. In such voids, the effect of cosmology would be 
dominated by the effect of peculiar velocities. If we found 
a way to characterize peculiar-velocity dominated voids, we 
could wisely exclude them from the analysis, as the signal- 
to-noise ratio for cosmological information with these voids 
would be low. In the next Subsection, we use the matched 
fraction of voids and the ellipticities to characterise the effect 
of peculiar velocities on cosmic voids. 


3.1 Matched fraction 

A void that is weakly affected by peculiar velocities will 
have similar properties in both simulations — a direct way 
to exclude the most affected voids is to consider the matched 
and unmatched voids. 

Taking the peculiar velocity catalogue as a basis for the 
comparison, unmatched voids from the non-peculiar sample 
are the most affected by peculiar velocities. For this analysis 
we considered ellipticity and radius of voids, to understand 
which properties of the unmatched voids are a feature of a 
velocity-dominated void. Figure 3 shows the matched and 
unmatched voids from the catalogues with and without pe¬ 
culiar velocities in the radius-ellipticity plane for both the 
HighRes and LowRes sample. The matching is worse for 
voids with radii lower than ~ 20 h“^Mpc for the HighRes 
sample and lower than ~ 35 h“^Mpc for the LowRes sam¬ 
ple, indicating that, when finding small voids, results can be 
strongly affected by peculiar velocities. This seems partic¬ 
ularly logical: we might have expected, a priori, that small 
voids are the most affected by changes in shape due to the ve¬ 
locities, which affects the way VIDE defines the Voronoi cells 
and selects them as belonging to a void. Additionally, since 
small voids are more likely to be found in higher-density en¬ 
vironments, they have higher ridges (Hamaus et al. 2014). 
Thus they are more affected by the peculiar velocities of 
the high-density structures forming the walls. The finding of 
these voids is affected by peculiar velocities, thus their fea¬ 
tures are peculiar-velocity dominated. Figure 3 shows that, 
while some properties of voids (such as ellipticity) are af¬ 
fected by peculiar velocities, the effect is not the same for 
all voids, and does not show any preferential feature that 
can be used to identify peculiar velocity-dominated voids. 
On the contrary, the radii of voids seem to allow to iden¬ 
tify unmatched voids among the smallest ones: properties of 
smaller voids are more likely to be dominated by peculiar 
velocities. 

When extracting cosmological information from the un¬ 
matched voids, the relative strength of peculiar velocities 
can lead to biases or systematics: the signal-to-noise is low 
due to the strong impact of velocities. With the aim of ex¬ 
tracting cosmological information from voids — for instance 
using Alcock-Paczyhski test — this consideration should be 




Figure 3. We represent in the top plot of each panel the matched 
and unmatched voids from the catalogues in the radius-ellipticity 
plane for both the HighRes and LowRes sample. Red crosses rep¬ 
resent unmatched voids; blue dots represent matched voids. The 
bottom plot of each panel shows the fraction of matched voids. 
Smaller voids can be significantly contaminated by peculiar ve¬ 
locities — they include a large fraction of the unmatched voids 
and are peculiar velocity-dominated (while here is no clear dis¬ 
tinction in ellipticity). We show that voids smaller than ~ 20 
h“^Mpc for the HighRes sample and smaller than 35 h“^Mpc 
in the LowRes sample are more affected by peculiar velocities. 
These limits in radius roughly correspond to twice the mean par¬ 
ticle separation. However, we also note that there are a population 
of well-matched, minimally-affected voids at all scales. 


taken into account, as it would greatly improve the result: 
low radii voids should be wisely excluded from the analysis 
to maximize the signal-to-noise. 

The bottom plot of each panel in Figure 3 shows the 
fraction of matched voids, confirming the guidelines we ob¬ 
tain from the top plot: for small voids the matching between 
voids without and with peculiar is less reliable. Thus when 
we find voids in redshift space, the smallest voids are the 
ones most affected by peculiar velocities and for which the 
matching is often worse, in that case the fraction of matched 
voids is lower. The fraction of matched voids is higher than 
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80% for voids bigger than 20 h“^Mpc for the HighRes sam¬ 
ple and 30 h“^Mpc for the LowRes sample. 

The use of cuts when working with voids can be a pow¬ 
erful tool to ensure a reduction of systematics errors. We 
point out that, depending on applications, the considered 
cut might vary. This does not apply to cuts aimed to avoid 
Poisson dominated voids (namely once the mean particle 
separation) but to other following cuts. 


3.2 Cuts and efficiency 

The previous Section showed that small voids are peculiar- 
velocity dominated objects: the use of a cut that excludes 
voids of small radius could reduce the effect of peculiar ve¬ 
locities on voids. For this purpose one would ideally need 
to find the trade-off between reducing intrinsic variance due 
to peculiar velocities and increasing statistical variance (due 
to the fact that we would remove voids). Nevertheless, the 
motivation for the cut might change, thus changing the cut 
itself. For example if we wanted to study the effect of pecu¬ 
liar velocities and eventually measure the peculiar velocity 
around voids, the voids with higher information would be 
the ones below the radius cut, i. e. small voids. 

From these considerations about unmatched voids, if 
the detailed structure of a void is necessary for the analysis, 
one should avoid using the smallest voids, to be more precise 
voids smaller than twice the mean particle separation^ since 
they may be contaminated by peculiar velocities. 

For the HighRes and LowRes samples, this would mean 
excluding voids smaller than 16 h“^Mpc and 30 h“^Mpc 
respectively. While this cut might seem drastic for current 
surveys (such as SDSS DR7) and might excessively reduce 
the number of voids, it can be a good prescription to be 
adopted for future applications of the Alcock-Paczyhski test. 

As Figure 4 shows, it is possible to further isolate un¬ 
matched voids by applying, additionally to the cut in radius, 
a cut based on the density contrast of voids (using the den¬ 
sity contrast as defined by ZOBOV and thus used by VIDE, 
which is given by the ratio of the minimum density along 
the ridge of the void versus the minimum density in the 
void). We see that unmatched voids (z. e. , those strongly 
affected by peculiar velocities) cluster in the small radii, low 
density contrast corner of the diagrams — which seems intu¬ 
itively consistent. Voids with low density contrast are shal¬ 
lower, thus the effect of peculiar velocities on such voids is 
stronger. 

As an example, we show in Table 1 and Table 2 the 
results of the example cuts shown in Figure 4: excluding 
voids with radii below twice the mean particle separation 
in radius and with density contrast below 1.15. With these 
cuts the remaining unmatched voids are ~5% for the high 
density sample and ~2% in the low density sample. We point 
out that the unmatched voids are legitimate voids, the only 
reason we consider their exclusion is because they are more 
affected by velocities. 

Nevertheless, for current datasets, these cuts also re¬ 
move many unaffected voids. In Table 1 and 2 we see that 
the example cut leaves as wished a very small number of 
peculiar-velocity dominated voids (remaining unmatched), 
but at the expense of generally reducing the number of voids 
and in particular excluding many of the good voids (the frac¬ 
tion of matched removal is high in both cases). 


3.0 

2.5 

P 

^^ 2.0 

1.5 
1.0 


10 ^ 10 ^ 

Mpc] 


3.0 

2.5 

, 0 ^ 2.0 

1.5 

1.0 


10^ 

[/i“^Mpc] 

Figure 4. We show with blue dots the voids matched between the 
peculiar velocity and no-peculiar velocity catalogues, and with red 
crosses the unmatched. In both the HighRes (left) and LowRes 
(right) cases we show with grey lines an example cut based on 
the radius and density contrast that removes a large fraction of 
unmatched voids (see Tables 1 and 2). The considered cuts only 
use information of the catalogue with peculiar velocities and are 
therefore applicable to data from real surveys. 




Cut 

%Match 

Removal 

%Unmatch 

Removal 

% Total 
Removed 

%Remaining 

Unmatch 

2Rmps 

20.0 

63.9 

28.6 

9.98 

l.lbdviDE 

39.0 

77.4 

46.6 

8.34 

Both 

44.6 

87.5 

53.1 

5.26 


Table 1: Example of removal proeess with euts for the high den¬ 
sity sample. 


To avoid this effect, we might want to consider cuts 
that avoid the exclusion of too many matched voids. A good 
trade-off between cutting too many voids and excluding as 
many non-matched voids as possible can be found using our 
matching algorithm: one can consider a cut minimising the 
number of removed matched voids and maximising the num¬ 
ber of removed unmatched voids. We perform a minimisa¬ 
tion taking into account these criteria in the parameter space 
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Cut 

%Match 

Removal 

%Unmatch 

Removal 

%Total 

Removed 

%Remaining 

Unmatch 

2-Rmps 

19.6 

67.0 

25.4 

5.39 

1.15Svide 

43.3 

86.3 

48.5 

3.24 

Both 

49.1 

92.5 

54.4 

2.02 


Table 2: Example of removal proeess with euts for the low density 
sample. 


of the density contrast and radius. More precisely, we define 
the efficiency ( 77 ) of the cut as the ratio of the percentage 
of unmatched voids removed to the percentage of 

the matched voids kept 

T\Tunmatch 

_ removed /r\ 

^ ~ ivfmatch 
^ kept 

Ideally, we would want the efficiency to be 1 (100% of un¬ 
matched voids removed and 100 % of matched voids kept), 
thus, to find the best parameters for the cuts in radius and 
density contrast, we minimize the function: 

f {Scut, Rent) = \l-r]\ (6) 

Figure 5 shows the exploration of different values for 
the cuts for the high density sample (top) and for the low 
density sample (bottom). 

We see an ideal zone in a curved region where the choice 
of parameters maximises the efficiency of cuts to master the 
effect of peculiar velocities on voids by excluding the most 
affected ones. In this zone the function / is minimized, thus 
maximizing the efficiency. 

Additionally, we show in Figure 6 {HighRes sample) 
and in Figure 7 {LowRes sample) the percentage of matched 
voids kept for each cut (top) and the percentage of voids 
dominated by peculiar velocities that are removed for each 
cut (bottom). Using values from the ideal zone shown in 
Figure 5 gives the best efficiency. A choice of parameters 
in this ideal zone typically gives for the HighRes sample 
~ 70% of unmatched voids removed at the expense of ~ 30% 
of matched voids removed; and for the LowRes sample a 
removal rate of ~ 75% for the unmatched voids and ~ 25% 
for the matched. For reference, in case the choice of values 
was not in that region. Figures 6 and 7 show the percentage 
of voids after any cut on radius or density contrast. 

Current and future surveys may particularly proht of 
the application of such cuts, to obtain catalogues of voids 
with reduced systematics, particularly adapted for the ex¬ 
traction of the cosmological signal. 


4 RELATIVE IMPACTS OF VELOCITIES ON 
VOID PROPERTIES 

4.1 Ellipticity 

For the population of matched voids, we examine the relative 
change in ellipticity due to velocities by comparing voids 
in the two catalogues (with and without velocities) for the 
HighRes and LowRes samples. The result of the comparison 
is shown in Figure 8 . 

We find an average reduction of the ellipticity of voids 
(change lower than 10% of the value) for both the HighRes 
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Figure 5. Maximisation of the efficiency (see equations 5 and 
6, we represent the function /) depending on the different values 
of cuts on density contrast and radius for the HighRes sample 
(top) and the LowRes sample (bottom). White shows the best 
efficiency zone for cuts. 


and LowRes samples, although there is significant scatter. 
For the HighRes Sample the ellipticity slightly increases for 
radii larger than ~ 40 h“^Mpc, these large voids are not 
flattened. Recently, Paz et al. (2013) used the redshift space 
distortion of the void-galaxy cross-correlation function to 
study the dynamical environment of voids. They found that 
smaller voids present an inner region in expansion, but their 
walls are in a collapse stage; while bigger voids are in expan¬ 
sion. Considering this scenario we should expect a flattening 
of small voids and an enlargement of large voids. 

The dynamical effects of the evolution of voids com¬ 
pete with the effect that peculiar velocities will have on how 
the void finder assigns cells to individual voids. Thus the 
percentage variation of ellipticity is a result of both these 
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Figure 6. For the HighRes sample: percentage of matched voids 
kept for each cut (top) and percentage of voids dominated by 
peculiar velocities (unmatched voids) removed for each cut (bot¬ 
tom). For reference, grey thick lines show the ideal area for cuts 
(see Figure 5): a choice of parameters in this ideal zone typically 
gives ~ 70% of unmatched voids removed at the expense of ~ 30% 
of matched voids removed (corresponding to ~ 70% of matched 
voids kept). 


Figure 7. For the LowRes sample: percentage of matched voids 
kept for each cut (top) and percentage of voids dominated by 
peculiar velocities (unmatched voids) removed for each cut (bot¬ 
tom). For reference, grey thick lines show the ideal area for cuts 
(see Figure 5): a choice of parameters in this ideal zone typically 
gives ~ 75% of unmatched voids removed at the expense of ~ 25% 
of matched voids removed (corresponding to ~ 75% of matched 
voids kept). 


processes, related to the dynamical evolution of voids and 
to the selection of cells. 

An interesting conclusion is then reached through this 
work: while voids are generally expanding and one might ex¬ 
pect peculiar velocities to elongate individual voids, consid¬ 
ering the full cosmic web and the process of finding voids we 
obtain a different result. Peculiar velocities lead to a thick¬ 
ening of the walls separating voids, thus the void is flattened 
between growing walls. 

The analysis comparing the catalogues with and with¬ 
out peculiar velocities gives us the amount of void ellipticity 
due to peculiar velocities. The percentage variation is cal¬ 
culated with respect to the ellipticity in redshift space: this 
point is fundamental, since it allows a direct correction to 
be applied to measured voids. It means that, if we take voids 
in redshift space and measure their ellipticity, less than 10% 
on average of this ellipticity is due to peculiar velocities — a 
correction of the ellipticity taking into account this value en¬ 
hances the robustness of the Alcock-Paczyhski test. Indeed 
recent applications of the test show that applying a correc¬ 
tion for the ellipticity (to take into account the observed 
constant systematic ellipticity reduction due to peculiar ve¬ 


locities) allows to extract the signal from real data (Sutter 
et al. 2014c). 


4.2 Centers and sizes 

We compute the relative radius between the peculiar velocity 
and no-peculiar velocity sample (see Figure 9). Interestingly, 
for voids bigger than twice the mean particle separation, 
the radius remains stable both in the HighRes sample and 
the LowRes sample, and for small voids the radii increase 
by no more than twenty percent. Thus, despite the change 
in shapes, as discussed above, voids preserve their average 
volume under distortions from peculiar velocities. 

We compare the relative distance of macrocenters, in 
Figure 10, defined as where d is the distance be¬ 

tween the macrocenters and Res is the effective radius of 
the void with peculiar velocities. We see that voids are dis¬ 
placed due to the effect of peculiar velocities. Once again the 
displacement is larger for small voids, but reduces to ~ 20% 
for larger voids. However, the change is in the line-of-sight 
direction, and so uses of the void macrocenter for ISW or 
lensing analyses are not much affected. 
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Figure 8. Percentage ellipticity change due to peculiar velocities 
for the HighRes sample (top) and for the LowRes sample (bot¬ 
tom) for matched voids. Peculiar velocities contribute on average 
to an increase in the sphericity of voids. We find an average reduc¬ 
tion of the ellipticity of voids at all scales (percentage ellipticity 
change lower than 10% of the value, grey dashed line), except for 
the largest in the HighRes sample, which are extended by a small 
percentage. The shaded regions are the Icr Poisson uncertainties. 


4.3 Radial density profiles 

We showed in the previous Section that individual small 
voids can be more affected by peculiar velocities. The stack¬ 
ing of voids can alleviate the effect: even if the shape of the 
void is not perfectly stable, the averaged density profile can 
average out fluctuations induced by random perturbations 
of the tracer particles. This has been suggested in the case 
of dark matter particle simulations by Lavaux & Wandelt 
(2012). We test this claim in a more realistic case using the 
HOD galaxy samples by comparing the density profiles of 
stacks in both cases (with and without velocity) for stacks of 
various radii. Using the technique described in Sutter et al. 
(2014b), we consider the co-centered density profiles: we saw 
that the void macrocenter is slightly displaced due to the ef¬ 
fect of velocities. Thus, we shift the centres of voids in the 
velocity sample to the macrocenter of each matched void 
and we build the radial profile around the void. 

The profiles presented in Figure 11 show that the den¬ 
sity profiles of stacks are not strongly affected by peculiar 
velocities, except for the smallest voids. We point out that 
the use of the co-centering technique is adapted here, since 




R^(i[h ^Mpc] 

Figure 9. Relative radius between the peculiar velocity and no- 
peculiar velocity sample. The shaded regions are the Icr Poisson 
uncertainties. For voids bigger than twice the mean particle sepa¬ 
ration the radius remains stable both in the HighRes sample and 
in the LowRes sample. 

we saw that there is a slight displacement of voids, but the 
shape and the radii of voids do not change significantly. We 
find similar results also without the use of the co-centering 
technique. When using real data, we can just consider the 
density profiles, without any need to re-center them, because 
the density profile is not changed but simply displaced. We 
also notice the slightly enhanced compensation wall for small 
voids, as described in previous analysis (Sutter et al. 2012b; 
Hamaus et al. 2014). 

As a conclusion, current cosmological constraints rely¬ 
ing on density profiles of stacks are only very mildly affected 
by peculiar velocities; as it is shown in Sutter et al. (2014c), 
where a constant offset can be used as a first approximate 
way to take into account their effect. Indeed the profiles are 
spherically-averaged, so line-of-sight distortions have a re¬ 
duced impact for the density profile. 

4.4 Abundances 

We compare the abundance of voids in mocks with and 
without peculiar velocities. We already saw that some voids 
present in the peculiar velocity catalogue cannot be matched 
with voids in the catalogue without velocities. Considering 
only this effect, the number of voids should be higher in the 
velocity case. But the measure of the void abundance is also 
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Figure 10. Relative distance of the position of the macrocenter 
of voids when adding peculiar velocities. The shaded regions are 
the Icr Poisson uncertainties. 


sensitive to the voids that are only in the non-peculiar veloc¬ 
ity sample but have been disrupted by peculiar velocities — 
which means that they are not present in the peculiar veloc¬ 
ity catalogue. For most observational probes these voids can 
be neglected since we will never be able to measure them 
(absent reconstructions, our measures are always in redshift 
space). 

Nevertheless for cosmological applications relying on 
the comparison of theoretical predictions for void abun¬ 
dances with observations, this effect of peculiar velocities 
also needs to be considered, since any distortions of the 
abundance functions due to peculiar velocities must be 
marginalized over to generate constraints on cosmological 
parameters. We show in Figure 12 the number of voids for 
the HighRes and LowRes samples; with and without pecu¬ 
liar velocities. We observe that peculiar velocities impact 
the smallest voids: in the HighRes sample, the presence of 
velocities results in a reduction of the abundance of voids. 
This can be explained by the disruption of small voids by 
velocities. In the LowRes sample this effect is still present, 
but at a minor level, because with a low density sample the 
sensitivity to small voids is lower. 

In the next Section we comment the application of the 
considered cuts to real data, and we conclude with general 
guidelines to master the effect of peculiar velocities in void 
applications. 


5 CONCLUSIONS 

With the increase of applications using cosmic voids, the 
modelling of systematic effects affecting voids becomes cru¬ 
cial. The main systematics one has to deal with when ex¬ 
tracting cosmological information from voids is the presence 
of the peculiar velocities of galaxies. The analysis presented 
in this paper allows one to tackle the problem of peculiar ve¬ 
locities by studying their effect on void properties through 
a detailed comparison of two void catalogues obtained from 
mock galaxy samples with and without the presence of pe¬ 
culiar velocities. 

Overall, voids are only mildly affected by peculiar veloc¬ 
ities. For most applications that use stacked voids density 
profiles, abundances, macrocenters, or effective radii (and 
thus volumes), the effects are negligible. In contrast, indi¬ 
vidual void ellipticities suffer a consistent ~ 10% reduction 
(except for the largest voids which suffer an extension). 

Additionally, peculiar velocities can strongly disrupt a 
population of small voids. Thus care must be taken when 
attempting to include them in an analysis. Voids with ra¬ 
dius below twice the mean particle separation and with low 
density contrast are more likely to be contaminated. Since 
a population of well-matched voids persists at all scales and 
density contrast, a trade-off must be considered for cuts 
based on these properties. Aiming to preserve an accept¬ 
able number of voids for the extraction of the signal, it is 
possible to use the best values determined in Section 3.1. 

We stress that it is certainly possible to include all iden¬ 
tified voids in a cosmological analysis, but the effects of pe¬ 
culiar velocities must be accounted for in the subsequent 
calculations. If instead we wish to generate a catalogue of 
voids with minimal systematic uncertainties, the proposed 
cuts can be directly applied to voids from real data. For ex¬ 
ample, applications of the Alcock-Paczyhski test to stacked 
voids from current data might profit from the exclusion 
of the smallest voids (for which the cosmological signal is 
strongly washed out by velocities): at the expense of the sta¬ 
tistical weight from small voids, there is a gain in reduced 
shape variance when limiting the analysis to only larger 
high-contrast voids. The cutting methodology discussed in 
this paper can already be applied within the VIDE toolkit. 

We have analysed the effect of peculiar velocities on 
void properties (radius, ellipticity, density contrast, density 
profiles and abundances) which affect the way we measure 
voids. We have shown that voids are indeed unique places in 
the cosmic web: they are only slightly affected by peculiar 
velocity distortions, enabling a cleaner examination of any 
underlying cosmological signal. 

However, in order to achieve sub-percent level preci¬ 
sion measurements of cosmological parameters, residual sys¬ 
tematics must be taken into account. We have shown that 
the systematic contribution from peculiar velocities is small 
(and, in many cases, negligible): we considered the applica¬ 
tion of optimal cuts to reduce the impact of velocities for the 
extraction of cosmological information from voids in current 
and future surveys. The work presented in this paper is a 
first step towards the understanding and modelling of sys¬ 
tematics in a fully self-consistent treatment of cosmic voids 
as cosmological probes. 
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Figure 11. The figure shows the density profiles of stacks of different sizes for voids with and without peculiar velocities. For both 
samples, from left to right the stacks are for 15-20, 25-30, 35-40, 45-50 h“^Mpc radii. Top row is the HighRes sample: the density profile 
is less affected by peculiar velocities while increasing in radius. Bottom row shows the LowRes sample: we observe the same trend as for 
the HighRes sample, but the effect of peculiar velocities is washed out at larger radii. 



[/i-i Mpc] 



Figure 12. Void abundances for the HighRes and LowRes sam¬ 
ples, with and without peculiar velocities. 
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